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In  a  bistatlc  SAR  system,  iso-range  and  iso-Doppler  contours  are 
orthogonal  only  for  zero  bistatic  angle.  The  range  gradient  diminishes 
as  the  bistatic  angle  is  increased.  When  polar  format  storage  is  used, 
these  effects  result  in  increasing  smearing  (defocusing)  and  geometric 
distortion  of  the  imagery  as  the  bistatic  angle  is  increased.  Hence, 
such  formatting  is  not  matched  to  two-dimensional  transform  processing. 

A  solution  to  the  bistatic  SAR  image  processing  problem  is  provided  by 
a  new  technique  called  coherent  unskewing  and  reconstitution  logic  (Curl 
The  matched  processing  performance  of  the  Curl  technique  is  demonstrated 
using  computer  simulation  of  received  data  from  a  bistatic  SAR  system 
with  two-dimensional  geometry.  The  resulting  imagery  is  sharply  focused 
and  is  totally  free  of  geometric  distortion. 

Tx 

With  three-dimensional  system  geometry,  the  use  of  Curl  with  two-dimen¬ 
sional  data  storage  and  processing  (with  a  fixed  data  storage  plane 
corresponding  to  the  synthetic  aperture  center)  may  provide  a  projected 
image  with  acceptably  small  distortion.  Tradeoff  curves  are  given  for 
estimating  the  orientation  of  the  projection  plane  as  a  function  of 
system  geometry  and  other  pertinent  parameters. 
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EVALUATION 


Airborne  Instruments  Lab  (AIL)  was  successful  in  removing 
radar  image  distortion  for  bistatic  radar  configuration.  They 
demonstrated  a  technique  which  records  data  in  such  a  way  as  to 
remove  the  "doppler  walk"  through  resolution  cells  when 
generating  an  image  from  radar  data.  This  technique  will  have 
application  in  tactical  target  identification  when  the  target 
is  being  imaged  by  more  than  one  radar  and  will  reduce  the 
amount  of  aspect  angle  change,  or  time  on  target,  frequently 


recruired.  This  effort  is  related  to  TPO  R2E. 


recrair 

L 


DANIEL  *L .  TAURONEY 
Project  Engineer 


1.  INTRODUCTION 


Bistatic,  or,  more  generally,  multistatic  radar  has  certain  advantages  compared  to  mono¬ 
static  radar.  A  key  advantage  is  the  added  security  it  provides  to  the  receiving  platform  in  the 
modern  electronic  warfare  environment.  (By  eliminating  the  transmission  of  electromagnetic  energy, 
either  completely  or  at  least  in  the  direction  of  the  target  area,  its  vulnerability  to  antiradiation 
missiles  (ARM's)  is  essentially  eliminated.)  Another  advantage  is  improved  target  signature  quality 
because  of  the  multiple  aspect  angles  observed.  However,  the  various  advantages  are  bought  at  the 
cost  of  generally  increased  complexity  of  the  multistatic  radar  system. 

Fine  resolution  synthetic -aperture  radar  (SAR)  imagery  is  a  potential  means  of  obtaining  target 
classification  and  identification.  For  monostatic  SAR,  the  polar  format  method  of  processing  data 
from  a  spotlight  radar  avoids  the  problem  of  range  walk,  thereby  facilitating  the  obtaining  of  fine 
resolution  over  a  usefully  large  area.  For  bistatic  SAR,  however,  the  nonorthogonality  of  the  lines 
of  constant  range  and  constant  Doppler  does  not  permit  comparable  results  with  polar  format  process¬ 
ing,  unless  the  transmitter  and  receiver  are  nearly  collinear  from  the  target  area  (that  is,  when  the 
bistatic  angle  is  very  small).  A  means  for  achieving  the  appropriate  data  formatting  has  been  de¬ 
vised  by  AIL.  The  method  is  referred  to  as  Coherent  Unskewing  and  Reconstituting  Logic  (Curl). 

The  processing  technique  has  been  recently  demonstrated  under  an  1R&D  program,  using  optical 
data  recording  and  processing  techniques  (by  additions  and  modifications  to  AIL's  Mini-Spotlite  radar 
data  recording  and  processing  equipment).  The  objective  of  the  program  reported  herein  was  to  ex¬ 
tend  the  techniques  to  provide  digital  processing  solutions  both  for  the  case  where  the  transmitter  and 
receiver  are  in  the  same  plane  as  the  field  of  view  and  for  the  more  general  case  where  they  are  not 
in  the  plane  of  the  field  of  view. 

The  processing  solutions  for  the  two-dimensional  case  are  demonstrated  on  small  digital 
arrays  of  synthetically  generated  bistatic  radar  data,  using  a  synthetic  aircraft  target  consisting  of  a 
small  number  of  point  targets.  The  Curl  processing  demonstrates  the  formation  of  well-focused 
imagery,  exhibiting  fine  resolution,  low  side  lobes,  and  freedom  from  geometric  distortion.  The 
limitation  on  these  aspects  of  the  radar  imagery  have  been  found  to  be  equivalent  to  those  for  mono¬ 
static  radar  imagery  with  comparable  radar  and  geometric  parameters,  when  using  polar  format 
processing.  That  is,  insofar  as  processing  limitations  are  concerned,  if  sufficient  transmitted  band¬ 
width  and  angular  aperture  is  available,  Curl  processing  can  provide  imagery  that  is  fully  the  equiva¬ 
lent  of  that  provided  by  polar  format  processing  of  monostatic  radar  data.  However,  it  must  be 
recognized  that  there  are  important  other  differences  in  a  practical  radar  system.  In  particular, 
target  properties  (such  as  the  bistatic  radar  cross  section  compared  to  the  monostatic  radar  cross 
section  and  bistatic  shielding  compared  to  monostatic  shielding)  will  generally  result  in  imagery  dif¬ 
ferences  that  are  unrelated  to  the  processing  algorithms  used.  Further,  in  actual  SAR  systems,  the 
increased  transmitted  bandwidth  necessary  to  support  a  specified  range  resolution,  as  bistatic  angle 
is  increased,  may  be  impractical  or  unaffordable.  Similarly,  where  only  the  transmitter  or  receiver 
has  tangential  motion  with  respect  to  the  target,  such  tangential  motion  provides  a  synthetic  aperture 
only  half  that  in  the  corresponding  monostatic  case.  Thus,  while  Curl  processing  provides  an  elegant 
and  exact  solution  to  the  processing  of  synthetic  aperture  data  from  two-dimensional  bistatic  radar 
system  geometries,  it  is  not  a  panacea  for  the  other  difficulties  of  bistatic  synthetic  aperture  radar . 

With  respect  to  the  general  three-dimensional  bistatic  SAR  geometry,  where  the  departure 
from  the  two-dimensional  case  is  slight,  use  of  an  imaging  plane  that  represents  the  average  or  mid¬ 
position  over  the  synthetic  aperture  provides  a  convenient  practical  solution  with  little  distortion. 

Analysis  of  the  three-dimensional  bistatic  geometry  has  led  to  expressions  for  the  orientation 
of  the  image  projection  plane.  From  these  generalized  expressions  a  series  of  tradeoff  curves  have 
been  calculated.  The  curves  can  be  used  to  aid  in  the  selection  of  transmitter  and  receiver  site 
geometry  or,  in  a  multistatic  synthetic  aperture  radar  system,  to  select  the  most  advantageous  trans¬ 
mitter/receiver  combinations  for  a  given  engagement. 

The  report  concludes  with  a  discussion  of  additional  investigations  that  could  lead  to  determining 
image  distortion  and  resolution  limits  and  also  verifying  improved  means  for  reducing  image  distor¬ 
tion. 


2.  CURL  PROCESSING  FOR  TWO-DIMENSIONAL  CASE 


2. 1  GEOMETRY  AND  CURL  STORAGE  EQUATIONS 

A  detailed  discussion  and  derivation  of  the  Curl  processing  concept  is  contained  in  a  proprie¬ 
tary  AIL  Technical  Note,  "Curl  Image  Processing  for  Bistatic  Imaging  Radar, "  J-D739,  July  1978. 

2.2  COMPUTER  SIMULATION  OF  TWO-DIMENSIONAL  BISTATIC  SAR  OPERATION 
2.2.1  Operation  Simulated 

The  generalized  two-dimensional  bistatic  SAR  geometry  to  be  simulated  is  shown  in  Figure  1A. 
Basic  assumptions  are  made  that  are  analogous  to  those  used  in  the  polar  processing  of  monostatic 
SAR  data: 

•  The  geometry  is  planar,  with  the  image  plane  determined  by  the  transmitter, 
receiver  and  center  of  the  field  of  view. 

•  The  transmitter  and  receiver  are  both  sufficiently  remote  from  the  field  of 
view  that  their  direction  from  any  point  in  the  field  of  view  is  constant  at  any 
instant  of  time.  Thus,  the  differential  range  of  any  target  component  from 
the  reference  point  is  measured  along  the  line  from  the  transmitter  (or  re¬ 
ceiver)  passing  through  the  reference  point  to  the  perpendicular  projection 
of  the  target  point  on  that  line  (as  shown  in  Figure  1A). 

The  normal  operation  of  a  high  resolution  synthetic  aperture  radar  with  digital  data  recording 
and  processing  equipments,  when  imaging  a  target  area,  would  proceed  as  follows: 

•  Recording 

•  For  each  transmitted  pulse  the  received  coherent  video  data  would  be 
quadrature  sampled  at  a  fixed  rate. 

•  The  received  quadrature  video  data  samples  would  be  digitized  and 
recorded  on  magnetic  tape. 

•  The  set  of  samples  received  from  each  particular  transmitted  pulse 
would  be  identified  in  some  manner  for  later  use  in  data  processing. 

(Normally,  the  identification  would  be  by  time;  however,  target  posi¬ 
tion  or  other  quantities  relatable  to  the  data  storage  angle  could  also 
be  used. ) 

Image  Processing 

•  The  recorded  data  would  be  reformatted  (from  the  form  in  which  it  was 
recorded)  to  a  rectilinear  form  suitable  for  application  of  a  two-dimen¬ 
sional  Fourier  transform.  This  conversion  could  be  performed  in  one 
of  two  ways: 

•  For  each  point  in  the  rectilinear  data  array,  data  would  be  se¬ 
lected  from  the  nearest  corresponding  point  in  the  original 
data  array.  If  the  original  data  were  sufficiently  oversampled, 
the  resulting  errors  (of  data  at  the  nearest  neighbor  compared 
to  at  the  exact  point  desired)  would  be  negligible. 

•  For  each  point  in  the  rectilinear  data  array,  data  would  be  inter¬ 
polated  (normally  in  two  dimensions)  using  the  original  data  array. 

•  The  rectilinear  data  would  be  Fourier  transform  processed  in  the  two 
orthogonal  directions  (after  weighting,  if  desired,  for  reduction  of  side 
lobe  levels)  to  form  an  image  of  the  target  area. 
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B.  CURL  STORAGE  SURFACE 


Figure  1.  Generalized  Two-Dimensional  Bistatic  Geometry  and  Curl  Storage 
Surface 


2.2.2  Computer  Simulation  Description 


For  the  purpose  of  these  demonstrations  of  the  quality  of  imagery  obtainable  using  the  Curl 
bistatic  processing  process,  computer  simulation  of  the  received  video  signal  will  be  used.  This  is 
in  lieu  of  using  a  real  bistatic  radar  system  in  flight  or  on  a  test  range. 

Instead  of  assigning  some  particular  motion  of  the  transmitter,  receiver,  and  target  as  a  func¬ 
tion  of  time,  for  simulation  purposes  it  is  convenient  to  effectively  work  the  problem  backwards. 

That  is,  the  rectilinear  data  array  is  filled  (in  some  convenient  order)  by  determining  the  bistatic 
angle,  bisector  direction,  and  (fast)  time  corresponding  to  each  data  array  point.  Then  the  direction 
to  the  transmitter  and  receiver  are  determined  to  match  this  condition.  For  convenience  in  the  simu¬ 
lation  (and  to  correspond  to  the  condition  specified  in  the  Statement  of  Work)  the  receiver  (or  trans¬ 
mitter)  is  fixed  in  direction  (that  is,  has  no  tangential  velocity,  but  may  have  radial  velocity  with 
respect  to  the  target).  Then  the  moving  transmitter  (or  receiver)  direction  changes  at  twice  the  rate 
of  the  bistatic  angle  bisector  direction  change.  By  the  use  of  this  backward  approach,  only  the  de¬ 
sired  set  of  rectilinearly  formatted  data  is  obtained.  Thus,  the  problem  of  data  interpolation  and  its 
possible  errors  are  avoided  in  this  simulation. 

The  bistatic  SAR  input  data  (complex  coherent  video  data  samples)  are  simulated  by  the  follow¬ 
ing  steps  for  each  target  point: 

•  For  Og,  the  selected  transmitter  angle  at  the  center  of  the  synthetic  aperture 

(corresponding  to  the  center  of  the  array  in  which  the  simulated  data  is  stored), 
determine  the  X  and  Y  offset  components  of  the  center  of  the  simulated  data 
array  from  its  non-offset  position  [corresponding  to  the  selected  reference 
point  location  (XR,  YR)]  . 

•  Determine  0Q  for  target  point. 

•  Determine  0Q  for  initial  point  in  data  array  (X  =  1,  Y  =  1),  using  fixed  aQ  and 

•  Determine  T  for  (half  bistatic  angle). 

•  Determine  transmitter  and  receiver  range  component  from  reference  to  target 

point.  These  are  R  cos  (0q  +  aQ  -  0Q)  and  R  cos  (oq  -  0q),  respectively,  where 

R  is  the  radial  range  of  the  target  point  from  the  reference  point. 

•  Determine  range  frequency  and  phase. 

•  Determine  quadrature  components  of  signal  at  initial  data  point. 

•  Add  quadrature  components  of  any  other  target  at  initial  data  point. 

•  Determine  quadrature  components  of  signal  for  all  other  points  by  appropriately 
incrementing  both  X  and  Y  over  the  range  1  to  N  (where  N  =  number  of  data 
points  on  each  side  of  the  signal  data  array). 

In  order  to  scale  the  computer  simulated  data,  the  simulation  program  includes  the  following 
parameters  (fixed  for  a  given  version  of  the  program). 

•  Start  frequency  of  FM  waveform 

•  End  frequency  of  FM  waveform 

•  Duration  of  FM  waveform 

•  Resolution  (separation  of  image  array  points) 

For  a  given  simulation  run,  the  following  parameters  are  input  from  the  computer  terminal 
keyboard: 

•  Fixed  receiver  (or  transmitter)  angle 
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•  Transmitter  (or  receiver)  angle  at  center  of  the  synthetic  aperture 

•  Target  coordinates  (X,  Y) 

The  overall  base  program  for  simulating  bistatic  radar  data  and  applying  Curl  processing  per¬ 
forms  the  following  functions  in  the  order  given: 

•  Prints  program  name  and  characterization 

•  Loads  complex  data  array  with  zeros  to  ensure  no  carry  over  of  unwanted 
prior  data 

•  Accepts  keyboard  input  of  receiver  angle  (aQ),  in  degrees 

•  Accepts  keyboard  input  of  transmitter  center  angle  (Og),  in  degrees 

•  Accepts  keyboard  input  of  target  coordinates  (X,  Y),  in  pixels 

•  Loads  data  array  with  simulated  target  return  data  in  complex  form  by 
calculating  transmitter  angle  and  range  delay  that  correspond  to  each 
data  array  point  using  Curl  storage  equations 

•  Accepts  keyboard  input  for  additional  target  coordinates  and  adds  target 
data  to  data  array  until  no  further  targets  are  desired 

•  Optionally,  prints  input  data  array  (as  0  or  blank,  corresponding  to  <.0 
and  >  0,  respectively) 

•  Hamming  weights  and  FFT's  data  array  first  in  range  dimension,  then  in 
azimuth  dimension 

•  Prints  maximum  amplitude  (in  dB,  reference  1)  of  output  data  array 

•  Prints  quantized  amplitude  of  output  data  array  (coded  "#"  for  maximum 
amplitude  to  -0.3  dB,  coded  "9"  to  "0"  for  10  successive  3  dB  ranges, 
and  coded  "blank"  for  -30.3  dB  and  lower  amplitudes) 

A  number  of  modified  versions  of  this  program  have  been  created  for  producing  the  several 
related  plots  presented.  These  include  provisions  for: 

•  Loading  a  simulated  aircraft  target  that  is  stored  in  memory  as  eight  equal- 
amplitude  point  target  components  (in  place  of  manual  keyboard  entry  of  one 
or  more  target  points).  (Output  plots  use  only  the  symbols  "#",  "9",  ”8", 

"7",  and  "6"  to  avoid  image  confusion  by  including  the  lower  amplitude  levels. ) 

•  Processing  with  polar  processing  algorithm,  rather  than  Curl  algorithm. 

•  Processing  with  reduced  transmitting  bandwidth  fixed  at  a  multiple  (r:uch  as  2) 
of  that  necessary  to  obtain  the  desired  resolution  at  zero  bistatic  angle,  rather 
than  sufficiently  wide  to  support  the  desired  resolution  to  any  specified  bistatic 
angle. 

•  Simulating  fixed  bistatic  angle  operation,  rather  than  holding  either  receiver 
or  transmitter  angle  fixed  while  the  other  angle  is  varied  (producing  variable 
bistatic  angle  operation). 

2.2.3  Computer  Generated  Plot  Formats 

Two  types  of  computer -generated  illustrations  are  used  to  demonstrate  the  two-dimensional 
bistatic  data  processing.  These  are  described  in  turn. 

2.2.3. 1  Output  Data  Printout-This  is  a  computer -generated  plot  of  the  image  of  the  simulated  target 
obtained  by  a  two-dimensional  Fourier  transform  of  the  input  data,  A  (nearly)  square  array  of  32  x 
32  data  points  is  plotted;  12  levels  of  amplitude  are  displayed.  The  largest  magnitude  point  on  a  given 
plot  and  all  others  within  0.3  dB  of  that  magnitude  are  plotted  as  "#".  The  magnitudes  at  all  other 
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points  are  quantized  into  3-dB  bands  and  plotted  as  the  characters  "9"  through  "0"  and  "blank" 

(more  than  30. 3  dB  below  the  maximum  value).  To  facilitate  determining  the  coordinates  of  any 
points,  scales  (in  pixels)  surround  the  output  image  plot.  Note  that  the  orientation  of  the  data  array 
is  fixed  and  corresponds  to  viewing  the  scene  from  below  the  printout,  that  is  at  an  angle  of  90  de¬ 
grees  as  defined  in  Figure  1A.  With  proper  processing  the  image  should  be  invarient  with  any  trans¬ 
mitter  and  receiver  angles.  This  fixed  orientation  of  the  field  of  view  was  chosen  for  these  demon¬ 
strations  to  facilitate  detection  of  any  image  distortion  when  the  transmitter  and/or  receiver  aspect 
angle  is  changed. 

For  all  of  the  output  printouts  presented,  the  receiver  angle  is  0  degree  and  the  transmitter 
angles  are  varied  from  0  degree,  typically  in  45-degree  increments,  to  180  degrees.  Because  there 
is  mirror  symmetry  of  the  images  for  the  range  of  transmitter  angles  from  360  down  to  180  degrees 
to  those  of  in  the  range  of  0  to  180  degrees,  only  the  latter  range  is  presented.  Because  the  receiver 
is  fixed  at  0  degree,  the  bistatic  angle  is  equal  to  the  transmitter  angle  for  all  the  images  presented. 

2. 2. 3. 2  Input  Array  Printout -This  is  a  computer -generated  plot  of  the  data  array  used  to  explain 
the  phenomenon  of  range  resolution  being  dependent  on  bistatic  angle  for  finite  transmitted  bandwidth. 
The  plot  is  a  (nearly)  square  array  of  32  x  32  points.  These  points  are  coded  "blank"  to  represent 
normal  input  data;  they  are  coded  "0"  to  represent  the  absence  of  input  data.  The  latter  results  from 
the  range  scale  variation  applied  to  the  input  data  (by  the  Curl  algorithm)  to  compensate  for  the  reduc¬ 
tion  in  range  gradient  as  the  bistatic  angle  is  increased. 

2.2.4  Simulated  Imagery  for  Various  Processing  Algorithms  and  Motions 

2.2.4. 1  Curl  Processing — Extremely  Wide  Transmitted  Bandwidth-The  effectiveness  of  Curl  proces¬ 
sing  will  be  demonstrated  first  for  the  general  case  where  the  transmitter  and  receiver  effectively 
rotate  about  the  target  at  different  rates.  Without  loss  in  generality,  the  simulation  actually  assumes 
that  the  net  effective  rotation  rate  of  either  the  receiver  (as  suggested  in  the  Statement  of  Work)  or  the 
transmitter  is  zero.  As  a  further  simplification,  for  this  first  case,  it  is  assumed  that  the  trans¬ 
mitted  bandwidth  can  be  made  as  wide  as  is  necessary  to  support  the  desired  range  resolution  to  any 
specified  bistatic  angle  less  than  180  degrees. 

For  this  case  it  can  be  seen  from  Figure  2  that  for  all  combinations  of  transmitter  and  receiver 
angles  (and  the  resulting  bistatic  angles  selected)  the  resulting  imagery  is: 

•  Totally  free  of  geometric  distortion  (all  dimensions  are  faithfully  preserved). 

t  Sharply  focused,  with  fixed  resolution  determined  by  system  parameters, 
including  fixed  data  aperture  size,  independent  of  transmitter  and  receiver 
angles  (and  the  resulting  bistatic  angles) . 

It  might  be  expected  that  in  Figure  2G,  where  the  bistatic  angle  is  180  degrees  (at  the  center  of 
the  synthetic  aperture)  the  image  would  be  smeared  in  the  horizontal  (X)  direction.  This  would  be 
the  case  at  exactly  180  degrees,  as  the  range  gradient  is  then  zero.  However,  for  the  parameters 
chosen  for  this  series  of  Figures,  the  synthetic  aperture  is  about  3  degrees.  This  fact,  coupled  with 
the  extremely  large  bandwidth  used,  results  in  an  image  at  a  bistatic  angle  of  180  degrees  that  is 
nearly  identical  to  that  obtained  at  all  smaller  bistatic  angles. 

2. 2. 4. 2  Curl  Processing— Matched  Transmitted  Bandwidth-Next,  the  practical  application  of  Curl 
will  be  demonstrated.  In  general,  the  transmitted  bandwidth  cannot  be  made  as  wide  as  assumed  in 
the  first  case.  Instead,  a  narrower  transmitted  bandwidth  would  be  used,  adequate  to  support  a 
specific  resolution  at  some  bistatic  angle  that  is  appreciable  less  than  180  degrees.  Then,  using  a 
fixed  data  aperture,  constant  resolution,  of  the  specified  value,  is  obtained  at  all  smaller  bistatic 
angles.  Alternatively,  by  using  a  variable  data  aperture  that  is  matched  to  the  available  data  extent 
at  each  different  bistatic  angle,  the  system  resolution  becomes  finer  than  the  design  point  value  at 
smaller  bistatic  angles.  The  system  resolution  will  become  coarser  at  bistatic  angles  greater  than 
the  design  point  value.  Note  that  this  coarsening  of  the  resolution  is  not  a  defect  of  the  Curl  proces¬ 
sing;  it  is  the  inevitable  result  of  the  blstatic  geometry  itself.  This  variation  of  the  iso-range  con¬ 
tour  gradient  is  compensated  by  the  Curl  processing  algorithm,  thereby  eliminating  the  geometric 
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distortion  that  would  otherwise  occur  (as  is  shown  in  the  following  section  when  polar 
format  processing,  rather  than  Curl  processing,  is  used).  However,  this  also  reduces 
the  fast-time  extent  of  the  data  relative  to  the  data  array  (processing  aperture).  When  the 
data  extent  becomes  less  than  the  aperture  dimension  the  range  resolution  is  correspond¬ 
ingly  coarsened. 

To  demonstrate  this  situation,  the  transmitted  bandwidth  that  is  chosen  is  just  twice 
that  required  to  support  the  specified  resolution  under  the  most  favorable  condition  of  zero 
bistatic  angle.  Then,  at  bistatic  angles  greater  than  a  specific  value  (45  degrees  in  this 
case),  the  available  data  does  not  fully  fill  the  data  array  (processing  aperture)  and  the 
resolution  is  consequently  reduced. 

For  this  matched  bandwidth  case  it  can  be  seen  from  Figure  3  that,  for  all  combina¬ 
tions  of  transmitter  and  receiver  angles  selected,  the  resulting  imagery  is: 

•  Totally  free  of  geometric  distortion  and  properly  located  within  the 
field  of  view. 

•  Sharply  focused  with  resolution  becoming  less  fine  as  the  bistatic 
angle  increases  beyond  some  specific  value  (determined  by  the  chosen 
excess  bandwidth  factor). 

Note  thai  Figures  3A,  B  and  C  are  essentially  identical  to  the  corresponding  Fig¬ 
ures  2A,  B,  and  C  for  bistatic  angles  of  0,  45,  and  90  degrees,  respectively.  However, 
Figure  3D  (bistatic  angle  of  135  degrees)  is  noticeably  smeared  in  the  now-diagonal  range 
direction.  This  is  a  direct  result  of  the  greatly  reduced  range  gradient  at  135  degrees  com¬ 
pared  to  0  degree.  Similarly,  Figure  3E  (bistatic  angle  of  150  degrees)  is  further  smeared 
in  the  range  direction. 

The  reasons  for  the  reduced  range  resolution  as  the  bistatic  angle  is  increased  are 
shown  in  Figures  4  and  5.  Figure  4  shows  the  Curl  data  storage  surface  for  the  cases 
corresponding  to  Figure  3.  With  the  chosen  bandwidth  factor  of  2  the  processing  aperture 
is  filled  at  bistatic  angles  of  0  and  45  degrees.  However,  at  bistatic  angles  of  90,  135,  and 
150  degrees,  the  aperture  is  less  than  fully  filled.  These  cases  are  more  clearly  illustrated 
in  Figure  5,  which  shows  the  processing  aperture  (data  array)  for  those  cases  (0  =  90,  135, 
and  150  degrees)  illustrated  in  Figure  3  where  the  processing  aperture  is  less  than  com¬ 
pletely  filled.  Notice  that,  while  at  90  degrees  the  aperture  is  nearly  filled,  as  the  bistatic 
angle  is  increased  the  aperture  becomes  less  filled.  (At  90  degrees  38  of  the  total  1024 
data  points  in  the  array  are  unfilled;  with  data  at  135  degrees  654  are  unfilled;  while  at 
150  degrees  869  are  unfilled. 

If  a  specified  range  resolution  is  desired  at  a  specified  bistatic  angle  it  is  possible 
to  provide  it  by  appropriate  choice  of  the  bandwidth  factor .  Thus,  to  a  bistatic  angle  of 
almost  90  degrees,  a  bandwidth  factor  of  2  is  sufficient.  If  this  performance  is  desired  to 
a  bistatic  angle  of  135  degrees,  for  example,  then  a  bandwidth  factor  of  about  10  is  required. 
Figure  6  shows  the  desired  performance;  no  coarsening  of  resolution  to  135  degrees,  and 
slow  coarsening  of  resolution  at  bistatic  angles  greater  than  135  degrees.  Note,  however, 
that  at  145  degrees,  with  146  of  1024  data  points  unfilled,  the  resolution  coarsening  is  almost 
imperceptable.  At  155  degrees,  with  roughly  half  the  data  points  unfilled  (495  of  1024),  the 
resolution  is  about  half  the  initial  resolution. 

2. 2. 4. 3  Polar  Processing  of  Rotating  Target  (Fixed  Bistatic  Angle)  Bistatic  Data-For  the 
special  case  in  which  the  target  rotates,  or  equivalently  when  transmitter  and  receiver 
rotate  about  the  target  at  equal  rates,  the  bistatic  angle  remains  constant.  For  this  case, 


9 


polar  processing*  suffices  to  produce  an  Image  that  appears  undistorted  In  shape,  but  is 
reduced  in  linear  scale  and  displacement  from  the  reference  point.  The  image,  however, 
is  otherwise  undistorted  in  that  the  scale  change  is  the  same  in  die  orthogonal  directions; 
hence,  all  angles  are  preserved.  Thus,  if  scale  is  disregarded  or  corrected,  aircraft 
recognition  would  be  unimpeded.  Since  scale  would  appear  to  be  an  important  ingredient  in 
aircraft  recognition  algorithms,  it  would  appear  desirable  to  provide  scale  factor  correc¬ 
tion  to  the  image  or  in  the  processing.  (This  is  an  inherent  part  of' the  Curl  processing 
algorithm.) 

Figure  7  shows  polar  processed  rotating -target  imagery,  without  scale  factor  correc¬ 
tion.  Note  that,  in  general,  the  symmetrically  located  target  points  are  represented  by 
somewhat  different  clusters  of  numbers  (representing  amplitude  of  response).  This  is 
simply  the  result  of  the  coarse  quantization  used  and  the  fact  that  the  sampling  points  are 
not  symmetrically  located  on  the  amplitude  response  function.  Actually,  the  image  is  free 
of  any  distortion  except  for  the  scale  factor.  When  the  data  are  stored  prior  to  Fourier 
transform  processing  in  such  a  way  as  to  correct  the  scale  factor,  the  resulting  imagery  is 
identical  to  that  of  Figure  2  for  Curl  processing.  In  both  cases,  the  transmitted  bandwidth 
is  extremely  wide  so  that  no  significant  coarsening  of  the  resolution  is  apparent  as  the  bi¬ 
static  angle  is  increased. 

2 . 2 . 4 . 4  Polar  Processing  of  Variable  Bistatic  Angle  Data-The  foregoing  section  demon¬ 
strated  the  imagery  resulting  from  the  special  case  of  fixed  bistatic  angle.  For  this  case 
bistatic  polar  processing  produced  images  having  perfect  shape  but  incorrect  scale. 

Next  to  be  considered  is  the  general  case  in  which  the  bistatic  angle  is  varied  during 
the  synthetic  aperture.  Before  presenting  the  imagery  of  the  simulated  aircraft  used  to 
demonstrate  the  prior  cases  it  is  instructive  to  show  the  imaging  of  a  single  point  target. 

This  is  illustrated  in  Figure  8.  Figure  8A  is  the  locus  of  a  target  point  at  X  =  24  pixels  and 
Y  =  0  pixels.  In  the  upper  plot  the  locus  is  shown  in  a  single  32  x  32  pixel  field  of  view.  In 
the  lower  plot  the  locus  is  shown  "unfolded”  in  two  adjacent  tield-ot-view  areas  to  more 
clearly  illustrate  the  typical  form  of  the  locus— a  roughly  triangular  shape  with  the  sides  gen 
orally  curved  lines  that  form  cusps  at  the  vertices  of  the  "triangle. "  In  both  plots,  the  num¬ 
ber  adjacent  to  each  plotted  point  is  the  transmitter  angle,  which  is  shown  (generally  in 
15-degree  increments)  from  0  to  360  degrees.  For  other  values  of  X  the  locus  is  similarly 
shaped  and  scaled  so  that,  at  a  transmitter  angle  of  0  degree,  the  target  point  images  at  the 
correct  location  in  the  field  of  view. 

Similarly,  Figure  8B  is  for  a  target  point  at  X  =  0  pixels  and  Y  =  24  pixels.  (In  this 
case  the  third  side  is  a  straight  line  corresponding  to  the  transmitter  angle  of  180  degrees; 
in  Figure  8A  the  third  side  was  a  curved  line  and  180  degrees  was  represented  by  a  single 
point  on  the  locus.)  Also  correspondingly,  for  other  values  of  Y,  the  locus  is  similarly 
shaped  and  scaled  so  that  at  0  degree  the  target  images  at  the  correct  location. 

Figure  8C,  for  X=24  pixels  and  Y  =  24  pixels,  Illustrates  the  more  general  case  where 
neither  X  nor  Y  is  zero.  While  the  general  fcusped  triangle  shape  is  preserved,  the  locus  is 
no  longer  symmetrical  as  were  the  prior  cases  (about  the  axis  on  which  the  target  point  was 
located).  For  this  case,  if  the  X,  Y  coordinates  of  the  target  point  are  scaled  by  the  same 
factor,  the  locus  is  similarly  shaped  and  appropriately  scaled  so  that  it  passes  through  the 
correct  location  when  the  transmitter  angle  is  zero. 


♦Polar  processing  is  here  defined  by  angular  data  storage  equal  to  bistatlc  bisector  angle 
and  linear  range  storage,  as  in  monostatic  polar  processing. 
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The  significant  point  illustrated  by  Figure  8  is  that,  while  there  are  gross  similari¬ 
ties  in  the  shapes  of  the  loci  when  plotted  in  the  "unfolded"  format,  there  are  also  very 
significant  differences.  In  general,  the  loci  for  target  points  at  different  locations  within 
the  field  of  view  are  sufficiently  different  that  clearly  apparent  distortion  of  the  image  shape 
occurs  for  multipoint  targets.  This  is  illustrated  in  Figure  9  for  the  same  simulated  air¬ 
craft  target  used  previously.  Note  that  for  zero  bistatic  angle  no  distortion  is  apparent, 
while  at  45  degrees  bistatic  angle  the  distortion  is  quite  apparent.  At  larger  bistatic  angles 
the  distortion  is  sufficiently  severe  that  aircraft  recognition  is  severely  impeded  or  made 
totally  impractical. 

The  use  of  simple  polar  processing  (as  opposed  to  the  matched  processing  of  Curl) 
leads  to  a  further  difficulty  not  readily  apparent  in  Figure  9.  In  addition  to  distortion  of 
the  images,  displacement  of  the  image  within  the  field  of  view  also  occurs.  To  facilitate 
observation  of  the  image  distortion,  all  images  of  Figure  9,  except  that  of  Figure  9A  for 
zero  bistatic  angle,  have  been  offset  (by  interger  numbers  of  pixels)  in  both  the  X  and  Y 
directions  so  as  to  approximately  center  the  simulated  aircraft  image.  These  offsets  are 
shown  on  each  image  subtitle.  Thus,  it  is  seen  that,  in  general,  for  polar  processing  not 
only  is  the  image  distorted  in  shape  and  scale  but  it  is  improperly  located  within  the  field 
of  view  as  well,  further  compounding  the  problem  of  aircraft  recognition. 


Because  of  the  complex  form  of  the  image  distortion  and  position  error,  it  appears 
impractical  to  apply  a  transformation  that  would  correct  these  imaging  errors  after  the 
image  is  formed.  However,  as  was  demonstrated  in  Figures  2,  3,  and  6,  the  use  of  the 
Curl  algorithm  to  properly  store  the  data  prior  to  image  formation  easily  removes  the 
cause  of  the  image  distortion  and  position  error  and  results  in  properly  located  and  undis¬ 
torted  imagery. 
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Figure  6.  Simulated  Aircraft  Image— Curl  Processing— Bandwidth  Factor  =  10 
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Figure  8.  Locus  of  Single  Target  Points— Bistatic 
Polar  Processing 
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3.  BISTATIC  IMAGING  IN  THREE-DIMENSIONAL  RADAR  SPACE 


Two-dimensional  image  processing  of  the  two-dimensional  bistatic  SAR  geometry 
case  is  extended  in  this  section  to  the  three-dimensional  bistatic  SAR  geometry  case.  For 
two-dimensional  processing,  it  is  necessary  to  determine  the  performance  limitations  due 
to  the  radar  or  target  third  dimension.  The  primary  effect  of  the  radar  third  dimension  is 
in  determining  the  orientation  of  the  "image  projection-plane."  For  radar  imaging  of  air¬ 
craft  from  the  ground,  the  target  is  nearly  planar,  and  the  radar  designer  is  primarily  con¬ 
cerned  with  the  orientation  and  scale  of  the  image  relative  to  the  target;  that  is,  is  the . 
image  a  plan  view,  side  view,  or  front  view  of  the  aircraft  target?  The  projection  plane  is 
determined  by  a  combination  of  the  orientation  of  the  target  relative  to  a  reference  plane 
(containing  the  target  center,  transmitter  and  receiver),  and  the  orientation  of  the  projec¬ 
tion  plane  relative  to  the  reference  plane,  determined  by  the  target  (or  radar)  velocity  com¬ 
ponent  normal  to  the  reference  plane.  This  target  perspective  investigation  is  the  primary 
subject  of  the  remainder  of  this  report;  tradeoff  curves  are  presented  to  allow  convenient 
evaluation  of  the  effect  of  various  parameters  upon  perspective. 

If  the  target  is  not  planar,  or  if  the  target  plane  is  not  closely  parallel  with  the  image 
projection  plane,  then  slight  changes  in  the  image  plane  or  target  orientation  occurring  over 
the  interval  of  data  collection  (the  synthetic  aperture)  will  cause  loss  (or  limitation)  of  reso¬ 
lution.  The  resolution  loss  increases  as:  the  target  "height"  increases,  the  target  orienta¬ 
tion  change  increases,  or  the  image  plane  orientation  change  increases.  This  effect  is 
explored  in  the  Appendix  for  the  simple  case  of  a  target  undergoing  rotation  only .  The 
rotating  target  case  is  shown  to  be  highly  similar  to  that  for  monostatic  radar,  allowing 
easy  extrapolation  to  the  bistatic  radar  case. 

In  the  following  two  sections,  the  bistatic  Doppler  gradient  vector  (subsection  3.1) 
and  the  bistatic  image  projection  direction  (subsection  3.2)  are  defined.  A  derivation  of  the 
bistatic  Doppler  gradient  vector  and  bistatic  image  projection  direction  can  be  found  in  the 
Appendix.  In  subsection  3.3,  the  orientation  of  the  image  projection  plane  in  terms  of  eleva¬ 
tion  angle  is  defined.  Subsection  3.4  presents  design  tradeoff  curves  of  elevation  angle 
versus  three-dimensional  bistatic  geometry. 

3.1  DEFINITION  OF  THE  BISTATIC  DOPPLER  GRADIENT  VECTOR 

V(XwD> 

The  bistatic  Doppler  gradient  vector, - 0  -  is  determined  by  the  vector  sum  of 

— i 

of  the  receiver  Doppler  gradient  vector  y  ( *  R)  and  the  transmitter  Doppler  gradient 

vector,  V (w x x  X).*  The  V(^R  *  ft)  vector  is  determined  by  using  the  plane  containing 

the  receiver  field-of-view  (FOV)  center  and  relative  velocity  vector  between  the  FOV  center 
and  receiver;  it  is  identical  to  the  Doppler  gradient  vector  of  a  monostatic  SAR  at  the  re¬ 
ceiver  position.  The  transmitter  Doppler  gradient  vector  is  determined  in  the  same  way 
by  considering  a  monostatic  SAR  at  the  transmitter  location.  A  derivation  of  the  bistatic 
Doppler  gradient  vector  is  given  in  the  Appendix. 


*4iR  (O^  =  receiver  (transmitter)  rotation  rate 

R  (X)  =  receiver  (transmitter)  unit  vector  direction 


3.2  BISTATIC  IMAGE  PROJECTION  DIRECTION 

If,  for  the  bistatic  three-dimensional  geometry  case,  .the  usual  assumption  is  made 
that  the  target  ranges  to  the  receiver  and  transmitter  are  large,  then: 

*  Iso-range  contours  are  planes  orthogonal  to  the  bistatic  angle 
bisector  direction 

*  Iso-Doppler  contours  are  planes  orthogonal  to  the  bistatic 
Doppler  gradient  vector 

The  iso-range  and  iso-Doppler  contours  intersect  in  straight,  parallel  lines.  Any  point  on 
a  given  intersection  line  is  at  the  same  range  and  Doppler  as  any  other  point  on  the  same 
intersection  line.  Then,  all  points  on  a  given  intersection  line  will  image  into  a  single  point. 
The  intersection  lines  then  give  the  image  projection  direction  which  is  unique  in  a  given 
geometry.  Any  plane  orthogonal  to  the  projection  direction  is  definable  as  an  image  projec¬ 
tion  plane  and  the  projection  direction  is  then  the  normal  to  the  image  projection  plane. 

It  is  important  to  note  that  the  iso -range  contours  (planes)  are  not  orthogonal  to  the 
iso-Doppler  contours  (planes)  unless  the  bistatic  angle  is  zero  (monostatic  case). 

The  following  equations  (derived  in  the  Appendix)  give  the  direction  of  the  image 
projection  plane  normal,  IT,  for  the  indicated  SAR  geometries: 

*  Monostatic  (General  Target  Motion) 


N  =  M  x  Ot  x  M  (1) 

•  Bistatic  (Rotating  Target) 

N  =  B  x  «T  x  B  (2) 

•  Bistatic  (Translating  Target) 

N  =  (1  +  cos  0)  («R+«X)  -  (X  •  riR)  R  -  (R  •  X  (3) 

•  Bistatic  (Translation  and  Rotation) 

N  =  (1  +  cos  0)  (Or  +  Ox)  -  (' 0R  •  B)  k  -  (Ux  •  B)  X  (4) 

where, 

0  3  bistatic  angle 

_ *  ,  _ 

M  3  monostatic  SAR  position  vector  with  respect  to  target  (called  2X  in  Appendix  A) 
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B  =  bistatic  angle  bisector  direction 
=  target  rotation  vector 


X  =  transmitter  direction  (unit  vector) 

A 

R  =  receiver  direction  (unit  vector) 

(ttjj)  =  bistatic  rotation  vector  of  receiver  (transmitter) 


Comparison  of  the  general  monostatic  case  (equation  1)  and  the  bistatic  case  for  a 
rotating  object  (equation  2)  shows  that  the  normal  to  the  image  projection  plane  has  the  same 
form  with  the  monostatic  radar  vector  5!  replaced  by  the  bisector  vector  a.  For  a  bistatic 
radar  configuration,  with  constant  bistatic  angle  and  aspect  diversity  generated  by  target 
rotation,  a  monostatic  equivalent  of  the  bistatic  radar  can  be  derived  with  the  monostatic 
radar  located  on  the  bistatic  angle  bisector. 

Comparison  of  the  bistatic  translating  target  case  (equation  3)  or  the  bistatic  transla¬ 
tion  plus  rotation  case  (equation  4)  with  the  monostatic  case  (equation  1)  shows  that  the  image 
projection  plane  normal,  N,  for  the  bistatic  cases  is  different  in  form  from  the  monostatic 
case.  The  translating  (and  translating  plus  rotating)  bistatic  cases  cannot  be  simplified  to 
a  monostatic  model  in  which  the  bistatic  terminals  are  replaced  by  a  single  fixed  monostatic 
terminal. 


3.3  IMAGE  PROJECTION  PLANE  ELEVATION  ANGLE 


For  a  translating  object  or  bistatic  SAR  (refer  to  the  Appendix),  the  projection  plane 
elevation  angle,  0,  relative  to  a  reference  plane  (containing  the  target  FOV  center,  receiver 
and  transmitter)  is  given  by. 


tan  0  = 


flL 


A  -A 

X  -  o. 


R 


v<v 


(R  *  X)  cos  0/2 


(5) 


where, 


X(&)  =  transmitter  (receiver)  unit  direction  vector 

_ i  — * 

R  =  transmitter  (receiver)  rotation  vector  due  to  target  translation  only 

A 

B  =  bistatic  angle  bisector  unit  direction  vector 


0  =  bistatic  angle 


This  expression  allows  the  magnitude  and  direction  of  0  to  be  easily  separated  with: 


direction  =  B 


magnitude  =  0  =  tan 


-1 


a. 


X  - Ox •  R 


(0X  +  «R)  •  (R  *  X)  cos  0/2 


(6) 


The  magnitude  of  6  is  a  good  measure  of  the  three-dimensionality  of  a  bistatic  engagement 
and  is  presented  in  a  set  of  tradeoff  curves  in  subsection  3.4.  The  translating  target  case 
(no  target  rotation)  is  of  particular  interest  as  applicable  to  imaging  a  nonmaneuvering  air¬ 
craft  with  a  ground-based  bistatic  SAR. 

3.4  DESIGN  TRADEOFF  FOR  A  SPECIFIC  ENGAGEMENT  GEOMETRY 

The  specific  engagement  geometry  for  imaging  a  noncooperative  aircraft  with  a  ground- 
based  bistatic  synthetic  aperture  radar  will  now  be  considered  and  analyzed.  Geometric 
image  distortion  will  be  considered  in  terms  of  the  average  elevation  (tilt)  angle  of  the  image 
projection  plane  for  the  synthetic  aperture.  The  derivation  of  the  image  projection  plane  ele¬ 
vation  angle  given  in  this  section  is  equivalent  to  that  given  in  subsection  3.3,  but  lends  more 
readily  to  a  series  of  pictorials  to  give  insight  and  understanding  to  the  derivation. 

The  remainder  of  this  section  develops  tradeoff  curves  to  assist  the  system  designer 
in  placement  of  transmitted s)  and  receiver(s)  to  minimize  image  distortion.  Additionally, 
the  tradeoff  curves  can  be  used  in  an  operational  multistatic  system  to  choose  the  optimum 
transmitter  and  receiver  to  minimize  image  distortion. 

3.4.1  Coordinate  System 

Consider  a  ground-based  bistatic  SAR  with  the  transmitter  and  receiver  arbitrarily 
located  with  respect  to  a  FOV  containing  an  aircraft  to  be  imaged  in  flight.  Without  loss 
in  generality,  the  coordinate  system  can  be  oriented  so  that  the  receiver,  transmitter,  and 
FOV  center  fall  in  the  X-Y  plane.  In  addition,  the  transmitter  can  be  placed  along  the  X- 
axis  and  the  target  FOV  center  at  the  origin.  The  target  aircraft  velocity  vector  is  then 
generally  not  in  the  X-Y  plane  and  can  have  a  significant  component  of  velocity  in  the  Z-axis 
direction. 

Figure  10  gives  the  relative  positions  and  directions  of  the  receiver,  transmitter, 

FOV  center,  and  the  target  velocity  vector  in  this  coordinate  system.  Also  indicated  on  the 
figure  is  the  bistatic  angle,  0,  and  the  direction  of  the  bistatic  angle  bisector  (which  also 
falls  in  the  X-Y  plane).  A  general  rule  of  thumb  for  small  bistatic  angles  is  that,  to  form  a 
synthetic  aperture  for  a  given  engagement  and  image  resolution,  the  direction  of  the  angle 
bisector  must  change  direction  by  a  minimum  prescribed  angular  amount  with  respect  to  the 
target.  If  Figure  10  is  valid  at  the  arithmetic  average  position  of  the  bisector,  it  is  also 
valid  at  the  angular  center  of  the  synthetic  aperture.  Further,  if  the  coordinate  system 
origin  is  fixed  to  the  target  FOV  center  in  time,  the  transmitter  and  receiver  will  have  iden¬ 
tical  velocity  vectors  equal  to  the  target  velocity  vector.  The  transmitter  and  receiver  will 
then  be  in  the  X-Y  plane  only  at  the  center  of  the  synthetic  aperture  (unless  the  Z-com- 
ponent  of  velocity  is  zero,  in  which  case  the  transmitter,  receiver,  and  FOV  center  are 
always  coplanar).  The  X-Y  plane  will  hereafter  be  called  the  fiducial  plane. 
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3.4.2  Image  Projection  Plane  Derivation 


The  orientation  of  the  image  projection  plane  at  a  given  instant  of  time  is  determined 
by  the  bistatic  angle  bisector  direction  at  that  instant  of  time  and  by  the  resultant  trans¬ 
mitter/receiver  Doppler  gradient  vector  at  that  instant  of  time.  The  resultant  bistatic 
Doppler  gradient  vector  is  the  vector  sum  of  the  transmitter  and  receiver  Doppler  gradient 
vectors. 

3.4.2. 1  Resultant  Doppler  Gradient  Vector-The  receiver  Doppler  gradient  vector  can  be 
computed  as  if  the  receiver  were  a  monostatic  SAR.  The  direction  of  the  receiver  Doppler 
gradient  vector  is  then  simply  given  by  the  tangential  component  of  velocity  to  the  receiver/ 
target  direction  in  the  plane  formed  by  the  receiver  position  and  the  target  velocity  vector. 
The  plane  formed  by  the  receiver  position  and  the  target  velocity  vector  is  the  image  pro¬ 
jection  plane  for  a  monostatic  SAR  at  the  receiver  location  (hereafter  called  the  receiver 
monostatic  image  plane). 

The  transmitter  Doppler  gradient  vector  can  be  found  in  the  same  way  by  considering 
a  monostatic  SAR  at  the  transmitter  location  and  computing  the  tangential  component  of 
target  velocity  in  the  transmitter  monostatic  image  projection  plane. 

The  receiver  monostatic  Doppler  gradient  vector  is  formed  by  taking  the  tangential 
component  of  target  velocity  in  the  receiver  image  projection  plane  and  normalizing  by  the 
receiver  to  target  distance.  The  transmitter  Doppler  gradient  vector  is  computed  in  an 
identical  fashion  considering  the  transmitter  as  a  monostatic  SAR.  The  resultant  bistatic 
Doppler  gradient  vector  is  the  vector  sum  of  the  transmitter  rotational  vector  and  the  re¬ 
ceiver  rotational  vector. 

Figure  11  considers  the  transmitter  position,  target  position,  and  target  velocity 
vector  at  the  center  of^the  synthetic  aperture.  The  transmitter  monostatic  image  plane  is 
determined  by  V  and  X  and  is  indicated  on  the  figure  by  plane  ABCO.  The  locus  of  all  points 
perpendicular  to  the  transmitter  direction  at  the  origin  of  the  coordinate  system  is  given  by 
plane  DCEO  whiqh  in  this  case  happens  to  be  the  Y-Z  plane.  The  tangential  component  of 
target  velocity  V^,  in  the  transmitter  monostatic  image  projection  plane  is  the  projection  of 

V  on  the  intersection  of  planes ^ABCO  and  DCEO.  The  transmitter  Doppler  gradient  vector 
is  then  formed  by  normalizing  by  the  magnetude  of  X*  to  get  wTX,  *  as  indicated  on  the 

figure.  It  is  important  to  note  that  Figure  11  is  not  drawn  to  scale  with  respect  to  the  rela¬ 
tive  lengths  of  the  vectors  V  and  X  for  clarity  of  demonstration. 

Figure  12  considers  the  receiver  position,  target  position,  and  target  velocity  vector 
at  the  center  of  the  synthetic  aperture.  The  receiver  monostatic  image  plane  is  determined 
by  V  and  R  and  is  indicated  on  the  figure  by  plane  ABCO.  The  locus  of  all  points  perpendi¬ 
cular  to  the  receiver  direction  (which  is  in  the  X-Y  plane)  at  the  origin  of  the  coordinate 
system  is  given  by  plane  DAEO.  The  construction  line  R^,  which  is  perpendicular  to  the 
receiver  direction  in  the  X-Y  plang,  has  been  added  to  clarify  the  picture.  The  tangential 

component  of  the  target  velocity,  V_0,  in  the  receiver  monostatic  image  projection  plane 
__k  TK 

is  the  projection  of  V  on  the  intersection  of  planes  ABCO  and  DAEO.  The  receiver  Doppler 
gradient  vector  is  then  formed  by  normalizing  VTR  by  the  magnitude  of  R  to  get  «TR  as 

indicated  on  the  figure.  Again,  it  is  important  to  note  that  Figure  12  is  not  drawn  to  scale 
with  respect  to  the  relative  lengths  of  the  vectors  V  and  R  for  clarity  of  demonstration. 


_ T  _ \ 

♦In  Appendix  AO  is  used  instead  of 
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Figure  10.  Geometry  at  Center  of  Synthetic  Aperture 


X 
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Figure  11.  Construction  of  Vector 


Using  Figures  11  and  12,  the  transmitter  and  receiver  Doppler  gradient  vectors, 
and  1?TR,  are  determined  at  the  center  of  the  synthetic  aperture  to  be: 


0,  -JS-,  4r 

w  |X|  lx  I 
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TR 


v  -  R„  v.  1 

v  x  x  yJ 
‘  - 

IRT 


RxfRxVy-Ryvxl 

-  ,3 - 


|R|* 


IrI 


(7) 

(8) 


where, 

(Xx>  0,  0)  =  transmitter  position 
(R  ,  R  ,  0)  =  receiver  position 

*  y 

(0,  0,  0)  =  target  FOV  center 
(vx,  vy,  vz)  =  target  velocity 

|x|  =  magnitude  of  the  transmitter  direction  vector 
|R|  =  magnitude  of  the  receiver  direction  vector 


The  bistatic  resultant  Doppler  gradient  vector,  ui,  is  given  by: 


w  =  WTX  + 


03 


TR 


(9) 


Figure  13  gives  the  vector  sum  of  u?__  and  w_,Y.  The  z-component,  u  ,  of  w  will  be  of 

IK  1A  Z 

interest  later  and  is  indicated  on  the  figure. 

3. 4. 2. 2  Bistatic  Angle  Bisector  Direction -The  bistatic  angle  bisector  direction  at  the 
center  of  the  synthetic  aperture  can  be  determined  from  Figure  10  to  be  given  by: 


B  =  (cos  [0/2],  sin  [0/2],  0)  (10) 

where  -  bistatic  angle  formed  by  the  transmitter,  target  FOV  center  and  the  receiver.  B 
is  in  the  X-Y  plane  at  the  center  of  the  synthetic  aperture. 

3. 4. 2. 3  Image  Projection  Plane  Elevation  Angle -The  image  projection  plane  elevation 

angle  with  respect  to  the  fiducial  plane  at  the  center  of  the  synthetic  aperture  is  determined 

from  Figure  14.  From  the  figure,  the  resultant  Doppler  gradient  vector,  is  obtained  from 

Figure  1  along  with  its  z-component,  u  .  The  bisector  direction  vector,  T?,  is  obtained 

z 

from  Figure  1.  The  bistatic  image  projection  plane  is  then  determined  by  the  directions  of 
u  and  B,  and  is  indicated  in  Figure  14  by  plane  ABCO. 
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From  Figure  14, is  in  the  image  projection  plane  and  gives  the  distance  wis  from  Tf. 

The  vector  A  is  in  the  fiducial  plane  and  connects  the  intersection  of  w  with  the  intersection 

z 

of  D  and  Then  the  elevation  angle  of  the  image  projection  plane,  0,  as  indicated  on  the 
figure,  is  given  by: 


0  =  tan 


I A  | 


or  using  equations  7  through  10: 


0  =  tan"* 


_ vJlxl+lRl] 

cos  (/S/2)  [  |X|  (v  sin  0  -  v  cos  0)  -  |R|  v  ] 
x  y  y 


where, 

(v  vy,  vz) 

0 

|X| 

IrI 


target  velocity  vector 
bistatic  angle 

distance  transmitter  is  from  origin  (target  FOV  center) 
distance  receiver  is  from  origin  (target  FOV  center) 


Rearranging  terms  in  equation  11  gives: 


(ID 


0  =  tan 


-1 


Let, 


cos  (0/2) 


IRI  V 

_|X|  +  |R|  \  V2  z  ] 

|x|  +  |r|  vz 

(12) 


- -  =  &  then  •  • -  i-  =  1-Efor  0  <  E  <  1, 


|X|  +  |R|  IX I  +  [R 

and  equation  12  can  be  rewritten  as: 


0  =  tan 


-1 


cos  (0/2) 

(1  -  E) 

1  V  V  \  V  ■ 

—  sin  0  -  -2-  cos  0]  -e/ 

\  z  2  1  ZJ 

(13) 


*  Equation  11  is  shown  to  be  identical  to  equation  6  (subsection  3-3)  for  this  specific 
geometry  in  the  Appendix. 
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The  normalization  by  |X  |  +  |?t  I  has  the  advantage  of  reducing  the  relative  distances  the 
transmitter  and  receiver  are  from  the  FOV  center  to  a  single  scalar  R  which  is  bounded 
between  the  finite  limits  of  0  and  1.  R  (X)  are  measures  of  the  bistatic  imbalance  of  the 
system;  that  is,  the  difference  in  receiver  and  transmitter  distances  from  the  target. 

If  we  let 


then 


X  V 

—  =  v  cos  $  and  —  =  v  sin  $ 
v„  v„ 

z  z 


and  equation  13  can  be  rewritten  as: 


0  =  tan  * 


which  is  the  elevation  angle, 
plane. 


(_  v_cos  (/3/2)  [sin(/3-$)  -  R|sin(/3-#)  +  sin  $J]  J 

0,  of  the  image  projection  plane  with  respect  to  the  fiducial 


(14) 


From  equation  14  and  referencing  Figure  15,  v  is  the  ratio  of  the  target  velocity 
component  in  the  fiducial  plane  to  the  component  of  velocity  perpendicular  to  the  fiducial 
plane  and  is  a  measure  of  the  degree  of  deviation  from  the  two-dimensional  geometry  case. 

V2  2 

vx  +  Vy  becomes  more  significant  than  vz  and  the  geometry  ap¬ 
proaches  two-dimensional.  As  v  becomes  small,  v_  dominates  and  the  geometry  is  highly 

z 

three-dimensional . 


From  Figure  15  and  equation  14,  the  angle  4>  describes  the  direction  of  the  target 
velocity  vector  projection  on  the  fiducial  plane  with  respect  to  the  transmitter  direction. 
It  is  important  to  note  that  the  expression  for  0  could  be  renormalized  so  that  <t>  is  ref¬ 
erenced  to  either  the  receiver  direction  or  the  bistatic  angle  bisector  direction. 

Before  presenting  tradeoff  plots  for  receiver,  transmitter,  and  target  geometries, 
a  few  general  comments  are  in  order. 


* 

a 
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Figure  14.  Construction  of  Image  Projection  Plane  and  Elevation  Angle 


9-2490 


Figure  15.  Geometric  Parameter  Normalization 


The  image  projection  plane  elevation  angle  and  the  azimuth  angle  are  of  significance  in 
relationship  to  scatterers  on  the  target.  For  example,  consider  an  aircraft  flying  straight 
and  level  at  an  altitude  that  is  of  the  same  magnitude  as  the  distance  to  the  transmitter  or 
receiver.  The  principal  scatterers  of  the  aircraft  are  contained  in  a  volume  that  is  longer 
and  wider  than  it  is  high.  The  image  projection  plane  orientation  will  not  be  nearly  parallel 
to  the  target  principal  scatterer  plane  and  the  image  produced  will  not  be  a  top  (plan)  view 
of  the  aircraft.  This  may  hinder  pattern  matching  to  identify  the  aircraft.  But,  if  the 
average  orientation  of  the  image  projection  plane  is  known,  it  may  be  possible  to  transform 
the  stored  pattern  of  the  aircraft  to  be  an  equivalent  view  to  that  of  the  SAR  image. 

It  is  also  possible  that  a  choice  of  receiver  sites  may  exist.  If  the  orientation  of  the 
image  projection  p  lane  at  the  center  of  the  synthetic  aperture  is  computed  for  each  receiver 
site,  then  a  receiver  site  can  be  chosen  that  minimizes  the  mismatch  in  orientation  of  the 
image  projection  plane  with  respect  to  the  principal  target  scatterer  plane.  If  the  tilt  angle 
can  be  minimized,  then  the  bistatic  data  can  be  stored  using  two-dimensional  Curl  techniques 
with  storage  angles  referenced  to  the  geometry  at  the  center  of  the  synthetic  aperture. 

The  expression  for  image  projection  plane  elevation  angle  given  in  equation  14  is  not 
sufficient  to  totally  describe  the  orientation  of  the  image  projection  plane  with  respect  to  the 
fiducial  plane.  Figure  14  shows  that  the  azimuth  angle  of  the  image  plane  can  be  defined  as  the 
angle  between  the  transmitter  direction  and  the  line  formed  by  the  intersection  of  the  image 
plane  and  the  fiducial  plane.  From  the  figure  at  the  center  of  the  synthetic  aperture,  the 
azimuth  angle  is  one-half  the  bistatic  angle. 

It  is  also  important  to  note  that  the  elevation  angle  and  azimuth  angle  of  the  target  with 
respect  to  horizontal  can  easily  be  determined  from  the  target  altitude  (from  a  tracking 
radar)  and  can  be  added  vectorially  to  the  elevation  and  azimuth  angles  of  the  image  projec¬ 
tion  plane  with  respect  to  the  fiducial  plane  to  obtain  the  azimuth  and  elevation  angles  of  the 
image  plane  with  respect  to  horizontal. 

3.4.3  Tradeoff  Curves 


Form  the  previous  section,  equation  14  can  be  solved  for  R  (bistatic  imbalance)  to  give: 


~  tan  Q  cos  (ft/2) J  +  S*n 
sin  4>  +  sin  (0-4>) 


(15) 


where: 

R  =  ratio  of  receiver  range  to  sum  of  transmitter  plus  receiver  ranges 

v  =  ratio  of  velocity  component  in  the  fiducial  plane  to  velocity  component  perpen¬ 
dicular  to  fiducial  plane 

<t>  =  angle  between  the  transmitter  direction  and  the  velocity  component  in  the 
fiducial  plane 


=  bistatic  angle 

0  =  image  projection  plane  elevation  angle 


For  a  specified  maximum  image  projection  plane  elevation  angle,  0  ,  a  specified  v 

max 

and  a  specified  0,  R  can  be  plotted  versus  4  to  determine  values  of  relative  distances  that 
the  transmitter  and  receiver  are  from  the  target  for  which  the  image  plane  elevation  angle 
will  not  exceed  the  absolute  value  of  0 

max. 

Figures  16  to  18  are  plots  of  R  versus  4  for  the  indicated  values  of  bistatic  angle,  0, 
and  velocity  condition  v: 


Figure 

p  (degrees) 

x 

16A 

90 

5.0 

16  B 

90 

5.5 

16C 

90 

10.0 

16D 

90 

15.0 

16E 

90 

18.0 

16  F 

90 

30.0 

16G 

90 

— >  00 

17A 

45 

10.0 

17B 

45 

15.0 

18A 

135 

10.0 

18B 

135 

15.0 

Each  figure  contains  plots  that  bound  the  regions  of  R  and  4  in  which  0max  is  less  than  or 

equal  to  5,  7.5,  10,  and  20  degrees.  A  detailed  discussion  of  the  formats  of  the  curves  is 
in  order  prior  to  a  discussion  of  the  significance  of  the  curves.  To  this  end,  Figure  16E 
will  be  discussed  in  detail. 


From  Figure  16E,  the  regions  for  which  the  conditions  of  0  less  than  or  equal  to 

max 

5,  7.5,  10,  and  20  degrees  are  coded  as  follows: 

Horizontal  lines 


0  <5  degrees 

max  s 

0  „  <7. 5  degrees 
max  6 

0  _  <  10  degrees 
max  6 


0  <20  degrees 

max  6 


Vertical  lines 

Tilted  lines  with  positive  slope  viewing 
from  left  to  right 

Tilted  lines  with  negative  slope  viewing 
from  left  to  right 


The  region  for  which  0  £  20  degrees  also  includes  the  regions  for  0  <  io, 

max  max 

7.5,  and  5  degrees;  the  region  for  0_„  <10  degrees  also  includes  the  regions  for  0  < 

max  max 

7. 5  and  5  degrees,  and  so  on.  The  solid  unhatched  region  is  for  the  condition  where  ®max 
>  20  degrees. 
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The  transmitter  direction  is  indicated  on  Figure  16E  by  the  vertical  dashed  line 
drawn  at  ®  =  0  degree.  From  Figure  16  the  projection  of  v  on  the  fiducial  plane  is  in  the 
same  direction  as  the  transmitter  for  ♦  =  0  degree.  The  receiver  direction  is  indicated 
on  Figure  3-7E  and  is  along  the  direction  of  $  =  £  =  90  degrees. 

The  vertical  scale  of  Figure  3-7E  is  labeled  R  and  (1  -  X).  For  R.  >  0. 5,  the  receiver 
is  at  larger  range  than  the  transmitter.  For  R  <  0. 5,  the  transmitter  is  at  longer  range 
than  the  receiver.  An  example  is  used  to  demonstrate  the  usefulness  of  the  curves. 

Consider  an  imaged  aircraft  that  is  90  miles  from  the  transmitter  and  10  miles  from 
the  receiver.  R  is  then  computed  as: 


R  = 


10 

90  +  10 


0.1 


The  aircraft  vertical  component  of  velocity  is  30  knots,  its  horizontal  velocity  component 
is  540  knots,  and  its  heading  with  respect  to  the  transmitter  is  ®  =  -20  degrees,  v  is  then 
computed  as: 


v  = 


540 

30 


18.0 


The  bistatic  angle  is  0  =  90  degrees.  The  parameters  /S  =  90  degrees  and  v  =  5.0  prescribe 
the  curve,  Figure  16E,  be  used.  From  16E,  for  R  =  0. 1  and  #  =  -20  degrees,  the  image 
projection  tilt  will  be  about  5  degrees.  It  is  important  to  note  that  the  transmitter  and  re¬ 
ceiver  terminals  can  be  switched  in  locations  with  the  same  results. 

Examination  of  Figure  16E  with  respect  to  any  given  condition  for  9  shows  that 

max 

the  acceptable  regions  are  periodic  in  *  and  go  through  two  complete  cycles.  The  unhatched 

portion  of  the  curves  indicate  that  for  a  given  receiver/transmitter/target  geometry,  a  given 

value  of  9  is  not  necessarily  satisfiable. 
max 

In  viewing  the  R  versus  *  curves,  if  no  region  is  given  for  which  9  is  less  than 

max 

5  degrees,  for  example,  then  this  condition  cannot  be  met  for  that  particular  engagement 
geometry. 

From  Figure  16E,  for  tne  curve  marked  for  9  <  20  degrees  in  the  lower  right 

max 

corner  of  the  plot,  there  appears  to  be  no  right-hand  arrow  marker.  Since  4>=  +180  degrees 
is  identical  to  ♦=  -180  degrees,  the  arrow  is  completed  in  the  lower  left-hand  corner  of  the 
plot  at  the  same  value  of  R.  The  boundary  of  4>=  180  degrees  is  continuous  with  the  boundary 
of  #=  -180  degrees. 

It  is  important  to  remember  that  the  image  projection  plane  elevation  angle  is  mea¬ 
sured  with  respect  to  the  fiducial  plane.  The  elevation  angle  can  be  related  to  a  horizontal 
plane  (of  interest  for  straight  and  level  aircraft  flight),  if  the  height  of  the  target  is  known 
above  a  horizontal  plane  containing  the  transmitter  and  receiver  and  the  target  ranges  to 
the  transmitter  and  receiver  are  known. 
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A  survey  of  Figures  16A  through  16G  shows  the  progression  of  satisfactory  regions 

for  9  for  a  bistatic  angle  of  90  degrees  and  v  increasing  from  5  to  ®  (target  motion  in  the 
max 

fiducial  plane,  which  is  the  two-dimentional  case).  The  condition  of  9  <  5  degrees  does 

max 

not  give  any  appreciable  R/*  areas  on  the  curves  until  v  =  18.  At  v  =  5  (Figure  16A),  9 

max 

<  20  degrees  is  satisfied  for  only  a  small  percentage  of  the  plot.  With  v  doubled  to  a  value 
of  10,  considerably  more  conditions  of  R  and  $  satisfy  9max  £  20  degrees.  The  elevation 

angle  is  then  very  sensitive  to  v  or  the  component  of  target  velocity  in  the  fiducial  plane 
relative  to  the  component  normal  to  the  fiducial  plane. 

The  effect  of  the  bistatic  angle,  8,  can  be  examined  by  comparing  Figures  17A,  16C, 
and  18A  which  are  all  for  v  =  10  but  bistatic  angles  of  45,  90,  and  135  degrees,  respectively. 
A  comparison  of  the  figures  for  any  value  of  9max  shows  that  as  8  is  increased  the  area  for 

which  9max  can  be  satisfied  decreases.  A  similar  comparison  of  Figures  17B,  16C,  and 

18B  for  v  =  15  also  shows  a  decrease  in  area  with  increading  8-  In  general,  for  a  given  v, 
the  area  to  achieve  a  maximum  9  condition  will  decrease  with  increasing  8- 


An  important  application  of  R  versus  $  tradeoff  curves  is  to  aid  the  system  designer 
in  the  placement  of  the  transmitter(s)  and  receiver(s)  in  a  multistatic  system.  If  the  ex¬ 
pected  direction  of  target  aircraft  is  known  within  boundaries,  then  the  transmitter(s)  and 
receiver (s)  can  be  placed  to  minimize  the  image  distortion  due  to  image  projection  plane  tilt. 
In  addition,  given  target  acquisition  by  a  tracking  radar  and  a  choice  between  receiver  sites, 
the  receiver  site  can  be  chosen  to  minimize  image  plane  tilt  via  equation  14  as  well  as  to 
minimize  the  expected  length  of  time  to  collect  the  synthetic  aperture.  For  reasonable  eleva¬ 
tion  angles  (less  than  20  degrees),  two-dimensional  Curl  storage  referenced  to  the  center  of 
the  synthetic  aperture  can  produce  an  image  which  is  sufficiently  geometrically  correct  to 
identify  the  aircraft  through  pattern  recognition. 

3.5  CONCLUSIONS 

If  the  orientation  of  the  bistatic  image  projection  plane  can  be  held  close  in  angle  to 
the  principal  scatterer  plane  of  the  target  aircraft,  then  two-dimensional  Curl  storage/ 
processing  can  be  applied  to  produce  an  image  which  has  a  high  likelihood  of  being  recognized 
and  identified.  Curl  storage  can  easily  be  computed  from  the  radar  geometry  in  the  fiducial 
plane  defined  at  the  center  of  the  synthetic  aperture. 

The  orientation  of  the  image  projection  plane  has  been  expressed  as  a  series  of 
tradeoff  curves  which  can  be  used  to: 

•  Choose  transmitter  and  receiver  site  placement. 

•  Choose  a  transmitter/receiver  combination  for  a  given  engagement  in  the 
multistatic  case. 

The  elevating  angle,  9,  of  the  image  projection  plane  is  sensitive  to  the  following 
parameters  with  the  indicated  relationships: 

•  9  increases  with  increasing  bistatic  angle,  8 

•  9  increases  with  decreasing  parameter  v  (the  ratio  of  "hor^.ontal"  to 
"vertical"  component  of  target  velocity). 
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•  6  is  sensitive  to  aircraft  heading. 

The  Appendix  of  this  report  discusses  recommended  additional  investigations  to  determine 
image  distortion  and  resolution  limits  and  investigations  into  handling  those  cases  which 
two-dimensional  Curl  will  not  produce  an  image  adequate  to  be  identified. 


/ 
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APPENDIX  A.  GENERAL  ANALYSIS 


A.  1  FOUNDATIONS  OF  VECTOR  ANALYSIS 
A.  1.1  Two  Dimensions  (2-D) 

The  figure  below  shows  a  transmitter  at  0X,  receiver  at  i/> R,  and  target  at  a  direc¬ 
tion  (and  S  distance)  relative  to  the  field-of-view  (FOV)  center,  in  2-D: 


The  transmitter  and  receiver  phases  (0V  D)  relative  to  the  FOV  center  are: 

A,  x\ 

0X  3  ^  cos  (o  -  il>x)  0R  =  —■  cos  (o  -  tl>R) 

The  total  return  phase  is  0  =  0X  +  0R,  and  return  Doppler  is: 


wDsdt  ^  =  ^ 


Using 


j^cos  (o  -  =  sin  (o  -  0.)  ^ 


2irS 

“d3T 


4>X  sin  (a  -  0X)  +  i>R  sin  (o  -  t l>R) 


Let  a  -  and  0 .  =0..  Then, 


2irS 

WDS~ 


Ox  sin  9X  +  0R  sin  0R 


i 
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Drawing  Cl.  as  vectors  perpendicular  to  the  receiver  or  transmitter  directions  and 
in  the  plane: 


the  return  Doppler  is  computed  by  adding  the  components  of  Ox  and  C1R  in  the  direction 
of  the  target. 

Equivalently,  Ox  and  0R  can  be  vectorally*  added: 


0T 


where  the  superscript  T  denotes  that  the  Cl  vectors  are  drawn  in  the  unconventional  tangen¬ 
tial  direction  illustrated  previously.  In  Section  3  the  symbol  was  used  instead.  Then 

target  Doppler  is: 


WD"  X  S 


-T  1  -T 


S 


as  shown  in  the  next  illustration. 


*  Henceforth,  a  vector  is  denoted  by  a  bar:  (A),  its  length  by  the  letter  alone:  (A),  and  its 
direction  by  a  caret  (A),  the  "unit  vector. "  Scalars  have  no  embelishments:  (X). 
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aasaaaa 


Rotating  the  coordinates  shows  that  the  iso-Dop's  are  lines  perpendicular  to  0,  for  ex¬ 
ample:  nT  •  S  =  0T  •  S' 


A.  1.2  Three  Dimensions  (3-D) 

A  solution  can  be  obtained  in  three  dimensions  (3-D)  by  considering  the  receive  and 
transmit  processes  separately.  Repeating  the  discussion  in  paragraph  A.  1 . 1 ;  for  either 
the  transmitter  or  receiver  using  the  plane  containing  the  FOV  center,  target,  and  either 
the  transmitter  or  receiver;  then: 


as  before.  However,  the  transmitter  and  receiver  angles  6.  are  no  longer  in  the  same 

plane.  The  figure  and  discussion  on  page  54  still  applies;  however  it  must  be  remem¬ 
bered  that  Ov  is  in  the  plane  defined  by  the  target,  FOV  center  and  transmitter,  while 
•"  _ t 

0R  is  in  the  target,  FOV  center,  receiver  plane.  O  remains  the  3-D  vector  sum  of 

-T  -T 
°X  +  0R: 

2n  q  ,  //sT  =Tv  2ff  «  xT 
WD=  X  S  (°X  +  0R)=T  S  ° 

_T 

and  the  iso-Dop's  are  planes  perpendicular  to  D  .  Now  the  surfaces  of  constant  range 
are  planes  perpendicular  to  the  bistatic  bisector.  The  intersection  of  an  iso-Dop  and  an 
iso-range  surface  is  a  line  of  constant  range  and  Doppler.  The  image  projection  plane  is 
defined  by  (the  normal  is)  any  line  of  constant  range  and  Doppler.  Hence,  the  image  pro- 
-T 

jection  plane  contains  O  and  the  bisector. 

A.  1. 3  2-D  Doppler  Gradient 

Note  «T  is  the  Doppler  gradient:  V(S  •  (for  0^  a  constant  and  S  the  position 
vector,  customarily  called  r)  is  O  .  This  is  proven  as  follows.  Let  the  x-axis  be  aligned 
with  0T: 


z 


7u 


6u  A  bu 
t>x  x  +  by 


A 

z  by  definition 
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Therefore 


V 


=  V(s  •  oT)  =  ox  =  oT 


(A-l) 


A.  2  VECTOR  ANALYSIS 
A.  2.1  3-D  Dopple i  Gradient 

The  preceding  concepts  can  be  developed  generally,  using  vector  analysis.  *  Con¬ 
sider  a  rotating  rigid  body  (FOV)  and  stationary  radar. 


The  velocity  of  the  target  is  0  y  S.  The  component  directed  toward  the  receiver  or  trans¬ 
mitter  is: 


R  *  (0R  x  S)  or  X  •  (Ox  y  S) 


(A-2) 


The  Doppler  shift  is  the  sum  of  the  receiver  directed  component  with  the  transmitter 
directed  component: 


X« 


2n  *  R  *  (0R  y  S)  +  X  •  (Ox  y  S) 


(A- 3) 


The  bistatic  Doppler  gradient  is  the  vector  sum  of  the  monostatic  gradients: 


V[«  '  ®R 
V[r 


y  S)  +  X  •  (Ov  y  S) 


•V[i 


{0B  yS)  +V  X  -  (OvXS) 


(A-4) 


*  Henceforth  all  angular  velocity  vectors  are  drawn  conventionally. 


Since 


R  *  0R  x  S  =  nR  •  S  x  R  =  S  •  RxOp  (identity) 


S  •  R  X  Or 


+  V 


—  fi 


S  •  X  x  O, 


(A- 5) 

(A-6) 


which  is  the  vector  sum  of  the  two  monostatic  gradients.  Expanding  the  gradients  (as  was 
done  to  obtain  equation  A-l)  gives  the  Doppler  gradient: 


A.  2.2  Range  Gradient 

To  determine  the  range  gradient,  let  L  be  the  distance  from  the  transmitter  to  the 
target  to  the  receiver : 


L  »  |R  -  S|  +  |X  -  S|  (A-8) 

The  range  gradient  Vl  is: 

Vl=  V|R  -  S|  +  V|X  -  S|  (A-9) 

Since, 

|R-S|3R-S'R 

|X  -  SI  S  X  -  S  •  X  (A- 10) 

\7l  =  Vr-Vs  -  R  +  Vx-Vs*  X  (A-ll) 


Also,  since  distances  R  and  X  are  constant,  independent  of  target  position: 

Vr  =  o  Vx  =  o 
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So 


-  VlsVs  •  R  +  Vs  •  X  =  R  +  X  (A-12) 

If  the  bistatic  angle  bisector  is  B,  the  range  gradient  is: 

-Vl  =  B  =  R  +  X  =  (2  cos  |-)  B  (A-13) 

where  p  is  the  bistatic  angle  between  X  and  R. 

A.  2.3  Image  Plane 

Now  the  normal  to  the  image  projection  plane,  N,  is  the  vector  product  of  the  range 
and  Doppler  gradients:  * 


N  =  (R  +  X)  x 


Op  x  R  +  Ox  x  X 


(A-14) 


Note: 

The  image  projection  plane  normal  can  equivalently  be  obtained  from  the  rotation  vector  of 
the  bistatic  bisector.  The  equivalence  is  illustrated  below,  in  two  dimensions. 


dR  +  dX 
dt  dt 


The  rotation  of  the  bisector  (in  three  dimensions)  can  be  determined  as: 


B  v  ®  -  (ft  +  *) 


x 


dB 

dt 


but 


dB 

dt 


dR 

dt 


+  =  O 

+  dt  WR 


A 

R 


+  °X  y 


and 

6g  =  (R  +  X)  x  ^  Op  x  R  +  Ox  x  xj 
Then,  from  equation  A-14  N  ■  Og 
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Then,  from  equation  A- 14 


n  =  ob 

A.  2. 3.1  Rotating  Target  Case-Note  that  if  only  the  target  FOV  rotates.  0„  =  0..  hence- 

__  K  A 

forth  called  0T 

N  =  B  x  0T  y  B  (A-l  5) 


and 


-Vl  =  b 

-V(^a)-6T>'  6 


(A-l  6) 


In  the  monostatic  rotating  target  case  R  =  X 


-V 


-Vl  =  2X 

/Xujt 


2ir 


=  20t  y  X 


N  =  4X  y  Ot  y  X 


(A-l  7) 


which  are  identical  to  the  bistatic  gradients  and  projection  plane  normal  if  B  is  replaced 
by  2X.  Further,  the  Doppler  gradient  is  perpendicular  to  the  range  gradient  for  both  bi¬ 
static  and  monostatic  imaging  of  a  rotating-only  target.  Thus,  monostatic  polar  format 
processing  can  be  used,  and  identical  image  projection  will  be  obtained  even  if  the  radar 
is  bistatic.  These  results  verify  those  obtained  by  Walker1. 

For  a  rotating  target,  the  image  plane  tilt  is  best  calculated  relative  to  the  rotation 
axis.  From  equation  A-l  5: 


N  =  B  y  0T  x  B  =  (B  •  B)  0T  -  (B  •  0T)  B  (identity)  (A-18) 

_  A  ^ 

If  the  scale  factor  N  is  of  no  interest,  B,  0T  can  be  replaced  with  B,  0T 

N  =  6t-(B-Ot)B  (A-19) 


1  J.  L.  Walker,  "Range-Doppler  Imaging  of  Rotating  Objects,"  PhD  Thesis;  University 
of  Michigan,  1974. 
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-(8  •  nT)$ 


v 


The  angle  between  the  image  projection  plane  normal  and  the  rotation  axis  (e)  is 
90  degrees,  less  the  angle  between  0T  and  B.  Hence,  the  image  projection  plane  contains 

B  and  the  direction  perpendicular  to  B  and  O^  (that  is  B  x  0T).  The  projection  plane  is 

fixed,  and  resolution  loss  is  due  to  target  rotation  only. 


Now  the  rotation  error  is  O  sin  4.  The  target  locus  is  a  circle  of  radius  H  sin  e'.  These 
effects  are  discussed  by  Walker  for  monostatic  radar  and  his  analysis  can  be  applied  to  bi- 
static  radar.  The  following  conclusions  may  be  drawn: 

•  A  3-D  storage  system  eliminates  errors 

•  2-D  optimum  storage  requires  only  timing  the  receiver 
to  the  FOV  center  range  and  results  in  equal  error  for 
all  target  points  at  the  same  altitude 

•  Range  and  Doppler  resolution  are  both  diluted  equally  by 
(l/cos  e) 

Error  due  to  resolution  cell  (projection  plane)  motion  over  a  data  collection  interval  At  is 

p  -  H  sin  e  OAt 

For  bistatic  radar,  the  ultimate  resolution  limited  by  diffraction  of  the  synthetic  aperature 
is  additionally  diluted  by  the  bistatic  angle: 


P 


A 

20  At 


1 

*  cos  e' 


x 


1 

cos  0/2 


(A-20) 
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2  \ 

Minimum  resolution  is  obtained  when:  (0  At)  =  R  gin  ^  c08  ^  and  is  the  (A-21) 

greater  of 


I)  XH  tan  e'-1  _ X _ 

|  2  cos  ff/2  OAt  cos  e'  cos  0/2 


(A-22) 


The  simple  example  of  a  rotating  target  illuminates  the  "image  projection  plane" 
concept.  At  any  instant,  a  3-D  object  will  image  as  if  viewed  along  N.  If  the  object  is 
2-D,  parallax  js  meaningless,  and  the  image  will  be  foreshortened  but  otherwise  the  same 
regardless  of  N. 

For  larger  time  intervals,  the  processing  can  be  adjusted  to  keep  points  in  any  one 
plane  fixed.  The  locus  of  other  points  is  the  projection  along  N  to  that  plane. 

For  polar  processing  the  one  plane  is  defined  by  O.  Points  on  the  0  plane  image 
as  points,  points  off  the  ft  plane  as  circles  even  if  they  are  on  the  image  projection  plane. 
Only  the  direction  N  is  meaningful. 

A  key  question  is  "Why  is  the  one  plane  0?"  The  answer  is  that  if  processing  is 
perfect  for  the  one  plane,  it  must  be  perfect  for  any  target  angle  over  a  full  turn  of  0. 

If  the  one  plane  were  the  projection  plane,  a  particular  point  on  the  plane  at  one  instant 
would  be  off  it  at  the  next  instant,  hence  defocussed. 

A.  2. 3. 2  General  Case-Returning  to  the  more  general  case,  expanding  equation  A-14: 

N  =  X  y  (0X  x  X)  +  R  x  (0R  x  R)  +  X  x  (C2R  x  R)  +  R  x  (6X  y  X)  (A-23) 

Since  for  any  vectors 

U  X  (V  x  W)  -  (U  •  W)  V  -  (U  •  V)  W  (A-24) 

U  x  (V  x  U)  =  (U  •  U)  V  -  (U  •  V)  U  =  V  -  (U  •  V)  U  (A- 2 5) 
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the  cc  nponents  of  N  can  be  expanded: 

x  y  nx  x  x  =  ax  -  (nx  •  x)  x,  r  x  nR  x  r  =  nR  -  (nR  •  R)  R 

X  x  0R  x  R  =  (X  •  R)  nR  -  (X  •  0R)  R  =  COS  0  0R  -  (X  •  0R)  R 

r  x  6X  x  x  =  (x  •  R)  nx  -  (R  •  nx)  x  =  cos  ^  nx  -  (R  •  nx)  x 

Then: 


N  =  (1  +  cos  0)  (0R  +  nx) 


A  A 

0R  •  (R  +  X) 

_  J 

|  R  - 

nx  •  (R  +  x) 

xi 

i 

(A-26) 


(A-27) 


Equation  A-27  is  a  general  expression  for  the  direction  of  the  image  projection  plane, 

given  the  transmitter  (X)  and  receiver  (R)  directions,  and  the  apparent  target  rotations  as 

seen  from  the  transmitter  (0„)  and  receiver  (0D).  The  apparent  rotation  can  be  the  result 

X  K 

of  target  rotation  or  translation;  or  radar  translation;  or  any  combination. 

A.  2. 3. 3  Translating  Target  Case-In  a  desire  to  simplify  equation  A-27  further,  special 
cases  are  considered.  The  elementary  target-only  rotation  case  has  already  been  examined, 
and  shown  to  reduce  to  the  monostatic  rotating  target  case.  Now  the  case  of  translation  only; 
of  target,  receiver  or  transmitter,  or  any  combination  will  be  considered.  The  expression 
in  equation  A-27  or  A-14  can  be  used  directly,  however,  in  this  case,  a  simpler  expression 
can  be  obtained  by  calculating  the  orientation  of  the  projection  plane  relative  to  a  "tilted  ref¬ 
erence  plane"  defined  as  containing  the  target,  receiver  and  transmitter.  This  reference 
plane  has  greater  physical  significance  than  the  horizontal  (or  ground),  since  gravity  plays 
no  part  in  the  imaging  process.  This  choice  of  a  reference  plane  is  analogous  to  the  choice 
of  the  axis  of  rotation  as  a  reference  for  the  rotating  target  case. 

Now,  from  the  definition  of  a  scalar  product 

N  •  R  x  X  =  cos  e  sin  0  (A-28) 

where  e  is  the  image  plane  tilt  relative  to  the  reference  plane  R  x  X.  Applying  equation  A-28 
to  equation  A-27: 


N  sin  0  cos  e  =  N  •  (R  x  X)  =  (1  +  cos  0 )  (0..  +  0R)  •  (R  x  X)  - 


(B  •  0R)  R  •  (R  x  X)  -  (B  •  nx)  X  •  (R  x  X) 


(A-29) 


Using  an  identity, 

A  A  A 

R  •  R  x  X  = 


and  the  fact  that  the  vector  product  of  any  vector  with  itself  is  zero: 

fi  A  A 

X  -  R  x  R  =  0,  X  -  R  x  X  =  R  •  XxX=0 
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Then 


N  sin  8  cos  e  =  (flx  +  0R)  ■  (R  y  X)  (1  +  cos  /J)  (A-30) 

Likewise  from  the  definition  of  a  vector  product: 

|Ux  V|  =  UV  sin  y 

where  U  and  V  are  any  vectors,  and  y  is  the  angle  between  them.  The  axis  of  rotation 
of  y  is  the  direction  of  U  x  V.  Then: 

N  sin  8  sin  e  -  |N  x  (R  x  X)|  =  |(1  +  cos  0)  (0R  +  0^)  x  (R  x  X)  - 

(B  •  nR)  R  x  (R  x  X)  -  (B  •  Ox)  X  x  (R  x  X)  |  (A-31) 

Expanding  equation  A-31: 

N  sin  8  sin  e  -  |(1  +  cos  0)  0R  x  (R  x  X)  +  Ox  x  (R  x  X) 

(0R  •  B)  R  x  (R  x  X)  -  (Ox  •  B)  X  x  (R  x  X)  |  (A-32) 

Then  expanding  the  triple  vector  products  using  identity  (equation  A-24): 

N  sin  8  sin  e  =  (1  +  cos  0)  (0R  •  X)  R  -  (Op  •  R)  X  +  (Ox  •  X)  R  - 

(dx  •  R)  X  j  -  (0R  •  B)  (R  •  X)  R  -  (R  •  R)  X  - 

(0  •  B)  r  (X  •  X)  R  -  (X  •  R)  X  (A-33) 

XL  J 


Substituting 

^  A  A  A 

X  •  X  =  1  =  R  •  R 
and 

X  •  R  =  cos  0 

Then  collecting  terms,  equation  A-33  simplifies  to 

N  sin  8  sin  e  =  j  0R  •  X  -  Ox  •  R  +  (cos  8)  (Ox  •  X  -  0R  R)  (R  +  X) 


(A-34) 


(A-35) 
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Dividing  equation  A-36  by  equation  A-30  gives: 


nR  •  x  -  n  •  r  +  cos  $  (nv  •  x  -  nD  •  r) 

tan  e  =  — - -  -  - - s - .  * - s -  % 

(n  +  0R)  •  (R  y  X)  cos  /S/2 


(A-37) 


The  denominator  clearly  depends  only  upon  the  components  of  the  rotation  vectors  normal 
to  the  reference  plane.  The  numerator  in  fact  depends  only  upon  the  component  of  trans¬ 
mitter  rotation  perpendicular  to  the  transmitter  direction,  and  the  component  of  receiver 
rotation  perpendicular  to  the  receiver  direction,  which  can  be  verified  by  examining  equa¬ 
tion  A-14  from  which  equation  A-36  was  obtained.  Thus,  the  image  projection  plane  orien¬ 
tation  does  not  depend  upon  whether  or  not  the  target  is  rotating.  However,  if  the  target  is 
not  rotating: 


nx •  X  =  O  0R  •  R  =  O 


(A-38) 


and  equation  A-37  simplifies  to: 


nR  •  x - nx  *  H 

(Ox  +  Or)  *  (R  X  X)  cos  /S/2  B 

(A-39) 


as  the  image  plane  tilt  angle  e  relative  to  the  target/transmitter/receiver  reference  plane 
when  the  target,  transmitter  and  receiver  are  translating.  This  means  that  tilt  angle  mag¬ 
nitude  is  the  arctangent  of  the  scalar  fraction,  and  the  axis  of  tilt  is  the  bistatic  bisector 
direction  B.  This  result  is  explored  in  paragraph  3.  3  under  the  further  restriction  that 
only  the  target  translates,  or  equivalently  that  the  velocities  of  the  transmitter  and  re¬ 
ceiver  relative  to  the  target  are  identical. 

For  a  target  translating  at  velocity  v,  and  both  receiver  and  transmitter  stationary: 


R  y  v 


X  y  v 


(A-40) 


Then  equation  A-39  becomes 


R  x  v  •  X  X  y  v  ;  R 
R  "  X 


tan  e  =■ 


COSS  fev  .  (R  y  X) 


(A-41) 


tan  0  cos  0  = 


X  (R  y  v  •  X)-R(Xxv-R) 

X  (R  x  v)  •  (R  y  X)  +  R  (X  x  v)  •  (R  x  X) 


(A-42) 
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Since 


U-VxW  =  V-WxU  (identity) 


(A-43) 


tan  e  cos  p  = 


X  (v  •  X  x  R)  -  R  (v  •  R  x  X) 

X  (R  X  v)  •  (R  x  X)  +  R  (X  X  v)  •  (R  x  X) 


(A-44) 


As  in  Section  3,  defining  the  x,  y  plane  so  x  is  aligned  with  X  and  y  is  in  the  same  plane 
as  R  and  X 


tan  e  cos  p  = 


X  v  sin  p  +  R  v  sin  p 
z  z 


X  (R  X  v)  •  (R  x  X)  +  R  (X  x  V)  •  (R  x  X) 


(A-45) 


where  vx,  v^,  vz  are  the  components  of  v.  Since 


(T  x  U)  ■  (V  x  W)  =  (T  •  V)  (U  •  W)  -  (T  •  W)  (U  •  V)  (identity) 


(A-46) 


X  (R  x  v)  •  (R  x  X)  =  Xv  •  X  -  X  (R  •  X)  (v  •  R) 


=  Xvx  -  X  cos  P  (v  •  R) 


(A-47) 


R  (X  x  v)  ■  (R  x  X)  »  R  (R  •  X)  (v  *  X)  -  Rv  •  R 


=  R  cos  vx  -  Rv  •  R 


where 


v  •  R  =  v  sin  p  +  v  cos  p 

y  * 


Then  the  denominator  of  equation  A-44  becomes 


(A-48) 


(A-49) 


Xv  (1  -  cos  p)  -  Xv  sin  p  cos  8  -  Rv  sin  p  =  sin  p 


Xvy  cos  8  -  Rvy 


Xvx  sin  p 


(A- 50) 


Hence 


v  (X  +  R) 

tan  e  = - -  (A- 51) 

cos  p  [  x  j  vx  sin  p  -  vy  cos  0]  -  RvJ 

Equation  A- 51  is  identical  to  equation  11  (paragraph  3. 4.  2. 3)  which  was  derived  in  a 
slightly  different  fashion. 
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MISSION 

of 

Rome  Air  Development  Center 

RAVC  plans  and  executu  ruearch,  development,  tut  and 
selected  acquisition  pn.ogn.amt>  In  support  Of$  Command,  Control 
Communications  and  Intelligence  (C3 1)  activitiu.  Technical 
and  engineering  support  within  areas  of  technical  competence 
<s  provided  to  ESV  Program  Officu  (P0< s)  and  other  ESV 
elements.  The  principal  technical  mission  areas  are 
communications ,  electromagnetic  guidance  and  control,  sur¬ 
veillance  0)(  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


